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Experiaents on Metallic aaterlalo have shown that fatigue cracks reaaln 
closed during part of the load cycle under constant- and variable-aaplitude 
loading* These experinents have shown that crack closure is a significant 
factor in causing load-interaction effects (retardation and acceleration) 
on crack-gro%fth rates under variable-anplitude loading* 

The present paper is concerned with the development and application of an 
analytical model cyclic crack growth that includes the effects of crack 
closure. The rnodel was based on a concept like the Dugdale model, but was 
modified to leave plastically-defonaed material in the wake of the advancing 
crack tip. 

The model was used co correlate crack-growth rates under constant-amplitude 
loading, and to predict crack growth under aircraft spectrum loading on 
2219-T851 aluminum alloy sheet material. The predicted crack-growth lives 
agreed well with experimental data. The ratio of predicted-to-experimen*"^l 
lives r«mged from 0.66 to 1.48. These predictions were made using data 
from an ASTM E24.06.01 Round Robin. 

INTRODUCTION 

Experiments [l] on metallic materials have shown that fatigue cracks 
remain closed during part of the load cycle under constant- and variable- 
amplitude loading. The crack-closure concept has been used to correlate 
crack-growth rates under constant-amplitude loading [2] and is a significant 
factor in causing 1 Interaction effects on crack-growth rates (retardation 



and acceleration) under variable-anplitude loading. Crack-cloaure is caused 
by residual plastic defo mat ions remaining in the wake of an advancing crack. 
Measurements of crack-opening stresses are very difficult and have been made on 
only a few materials and for a limited number of loading variables. To develop 
the rationale for predicting crack growth under general cyclic loading, a 
mathematical model of crack closure must be developed and verified for constant- 
and variable-amplitude loading. 

The crack-closure phenomenon has been analyzed using two-dimensional, 
elastic-plastic, finite-element methods [3-6]. The finite-element analyses 
were shown to be quite accurate, but were very complicated and required large 
computing facilities. There have also been several attempts to develop simple 
analytical models of crack closure [3,7-llj. All of these models were based 
on a concept like the Dugdale model [12] or strip-yield model, but modified 
to leave plastically-deformed material in the wake of the advancing crack. 

Newman [3], Budiansky and Hutchinson [s], and Fuhring and Seeger [lO,ll] 
studied only the crack-closure behavior. But, Dill and Saff [?] and Hardrath, 
Newman, Elber and Poe [9] used the crack-opening stresses from the models to 
predict crack growth under spectrum loading. However, none of these studies 
considered the Influence of three-dimensional constraint on the crack-closure 
behavior. 

The purpose of the present paper is to develop and apply an analytical 
crack-closuio model that simulates plane-stress and plane-strain conditions. 

The present model was based on the Dugdale model, but was modified to leave 
plastically-deformed material along the crack surfaces as the crack advances. 
Plane-stress and plant-strain conditions were simulated by using a "constraint" 
factor on tensile yielding to account for three-dimensional effects. 
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The aodel was developed for a central crack in a finlte^^rldth speclaen 
that was subjected to uniformly applied stress. The nodal was essrcisad under 
constant-ai 4 >litude loading using various applied stress levels and stress 
ratios. The crack-opening stresses calculated fron the nodal were used to 
calculate Elber's effective stress-intensity factor range [l]. Experlaental 
crack-g*xNrth rates from 2219-T8S1 aluainun alloy sheet aaterial under constant- 
anplitvde loading were correlated with the effective stress-intensity factor 
range far a vide range of stress levels and stress ratios. The experlaental 
data were obtained from Chang and Stolpestad [l3]. An equation relating 
crack-growth rate to effective stress-intensity factor range, threshold 
stress- Intensity factor range, and fracture tou^ness has been developed for 
the complete range of crack-growth rates. The closure aodel was then used 
to predict crack growth in 2219-1851 aluminua alloy aaterial under aircraft 
spectrum loading [13]. The predictions trere made using the closure aodel 
under nearly plane-strain conditions. These predictions %fere made using 
data from an ASTM E2A.06.01 Round Robin. The purpose of the round robin was 
to compare various methods for predicting crack growth in 2219-T851 aluminua 
alloy center-crack tension specimens subjected to aircraft spectrum loading. 

LIST OF SYHBOLS 

dimensions for partially-loaded crack (k 1,2), m 
material crack-growth constants (k " 1,5) 
half-length of crack, m 
half-length of final crack, m 
half-length of initial crack, m 
half-length of starter notch, m 
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d 

E 

F 

K 



N 



n 



V 

W 


w 


1 


x.y 


X 


i 


a 


half-length of crack plus tensile plastic sons, ■ 

Young* s aodulus of elasticity, MPa 

boundary-correction factor on vi tress intensity 

1/2 

stress-intensity factor, NPa-a 

1/2 

stress-intensity factor at failure, MPa-n 

1/2 

■axlnun stress-intensity factor, MPa-n 

1/2 

elastic-plastic fracture toughness, MPa-n 

length of elenent (i) created by plastic deforaatlon, m 

fracture toughness paraoeter 

nunber of cycles 

nuaber of cycles predicted from analysis 

nuidl>er of cycles from test specimen 

total number of bar elements 

stress ratio (S /S ) 

' min max' 

applies stress, MPa 
maximum applied stress, MPa 
minimum applied stress, MPa 
crack-opening stress, MPa 
specimen thickness, m 
crack-surface displacement, m 
specimen width, m 

half-width of bar element at point 1, m 

Cartesian coordinates 

coordinate location for element i, m 

constraint factor, a ■ 1 for plane stress and a » 3 for plane 
strain 


A 



Ac* 

AK 


AK 


eff 


AK 


AK 


th 


AS 


eff 


n 



crack-growth increaaat over which is held constant 

1/2 

stress-intensity factor range* MPs-e 

1/2 

effective stress-intensity factor range, MPa-« 

1/2 

effective threshold stress-intensity factor range, MPar-e 

1/2 

threshold stress-intensity factor range, MPa-e 
effective stress range, MPa 

material constant, n ~ 0 for plane stress and n * v for plane 
strain 

Poisson's ratio 

length of tensile plastic zone, n 

stress on segment of crack surface, MPa 

flow stress (average between a and O ) , MPa 
® ys u 

yield stress (0.2 percent offset), MPa 
ultimate tensile strength, MPa 
length of compressive plastic zone, m 


ANALYTICAL CRACK-CLOSURE MODEL 

To calculate crack-closure and crack-opening stresses during crack 
propagation, the elastic-plastic solution for stresses and displacements in 
cracked body must be known. The crack-surface displacements, which are used 
to calculate contact (or closure) stresses during unloading, are influenced 
by plastic yielding at the crack tip and residual deformations left in the 
wake of the advancing crack. Upon reloading, the applied stress level at 
which the crack surfaces become fully open (no surface contact) is directly 
related to contact stresses. This stress is called the "crack-opening 
stress." 



Because there are no cloBed**fora solutions for elastle-^lastlc cracked 
bodies, siaple approxiaations »ist be used. As previously aentioned, the 
present model was based on the Dugdale model. In tlm proposed model, 
severs assumptions about the plastic sone, the material cyclic stress-strain 
behavior, the residual plastic deformations, and crack extension were made. 

Some of these assumptions are discussed in the following sections and the 
others, along with details, are discussed in appendix A. 

The model developed herein was for a central crack in a finite-width 
plate subjected to uniform applied stress, as shown in figure 1. The model 
was based on the Dugdale model, but was modified to leave plastically-deformed 
material in the wake of the advancing crack. The primary advantage in using 
this model is that the plastic-zone size and crack-surface displacements are 
obtained by superposition of two elastic problems. These two elastic problems, 
a crack in a finite-width plate subjected to either remote uniform stress, S, 
or to a uniform stress, a, applied over a segment of the crack surface, 
are shown in figure 2. The stress- intensity factor and crack-surface 
displacement equations for these loading conditions are given in appendix B. 

Figure 3 shows a schematic of the model at maximum and minimum apnlied 
stress. The model was compos J of three regions: (1) a linear elastic 

region containing a fictitious crack of half-length c + p, (2) a plastic 
region of length p, and (3) a residual plastic deformation region abong the 
crack surfaces. The physical crack is of half-length c. Regi^ n 1 was 
treated as an elastic continuum, and Che crack-surface displacements under 
various loading conditions are given in appendix B. Regions 2 and 3 were 
composed of rigid-perfectly plastic (constant stress) bar elements with a flow 
stress, which is the average between the yield stress, 
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ultlute tensile strength, The shaded regions in figure 3(a) and 3(b) 

indicate aaterial which is in a plastic state. At any applied stress level, 

the bar eleaents are either intact (in the plastic sons) or brcdcen (residual 

plastic deforeatlon) . The broken elenents carry coaqpressive loads only, and 

then only if they are in contact. The elenents in contact yield in conpression 

when the contact stress reaches -O . Those eleaents that are not in contact 

o 

do not affect the calculation of crack-surface displaceaents. To account for 

the effects of state-of-stress on plastlc-eone size, a constraint factor a 

was used to elevate the tensile flow stress for the intact eleaents in the 

plastic zone. The effective flow stress 0U7^ under simulated plane-stress 

conditions was o and under simulated plane-strain conditions was 30 . 

o o 

The constraint factor is a lower bound for plane stress and an approximate 
upper bound for plane strain. (Although the Dugdale model is not entirely 
suited for plane strain yielding, the size of the plastic zone at the crack 
ti ' and its influence on crack-surface displacements may be adequate.) These 
constraint factors were verified using elastic-plastic finite-element analyses 
of cracked bodies under plane-stress [6] and plane-strain conditions (by 
the author) . The procedure used to establish the constraint factor (a) is 
discussed later. 

The analytical crack-closure model, discussed in appendix A, was used 
to calculate crack-opening stresses, S^, as a function of crack length and 
load history. In turn, the crack-opening stress was used to calculate the 
effective stress-intensity factor range, as proposed by Elber, and, conse- 
quently, the crack-growth rates. 
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FATlGUE-CaMX GSOMTH BATE BQUATIOE 

Hm crack-growth equation proposed by BIber [2] etatea that the crack- 

b 

growth rate Is a power function of the effective stress-lnteneity factor range 
cmly-. Later» Hardrath» Newnan, Elber and Poe [9] ahoired that the power law 
was inadequate at hl^h growth rates approaching fracture. The results pre- 
sented herein show that it is also Inadequate at low growth rates approaching 
threshold. To account for these effects, the power lav was modified to 



and 


■ (S - S ) Ac F 
eff ' max o 


( 1 ) 


( 2 ) 

(3) 


(4) 


The crack-opening stresses, S^, were calcualted from the analytical closure 
model. Equation (1) gives the "sigmoidal" shape commonly observed when fatigue 
crack-growth rate data are plotted against stress- intensity factor range. In 
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tlM lnt«TMdlattt ranga of crack-growth rataa, aquation (1) la baalcally Ubar'a 

C9 

propoaad powar law, ^aff* conatanta to vara dataralaad 
to beat fit ai^arlaantal data for conatant-aaplltuda loading. 

The coafflclanta and vara dataralnad froa threabold data on 

the 2219-T8S1 alualnun alloy ahaat notarial fron reference [14]. The effective 

threabold atreaa-lntenalty factor range, waa determined froa the threabold 

atreaa-lntenalty factor range, as 


AK - U M . 
o tb 



(5) 




The coefficient is the elastic stress-intensity factor at failure 

1/2 

or cyclic fracture toughness. The coefficient was chosen to be 77 MPa-a 
1/2 

(70 ksi-in ) on the basis of the crack-growth tests in reference [13] 

(see appendix C) . 

The coefficients and C 2 were found from constant- amplitude rate 

data [ 13 ], after C^, and were determined, by using a least-squares 

regression analysis. The constant-amplitude correlations were made using 
values computed from the model with various constraint factors. 

It was found that an a of about 2.3 would give a good correlation. A 
summairy of the coefficients used to correlate the constant-amplitude data 
under nearly plane-strain conditions (a - 2.3 and n * 0) are as follows: 


- 1.764 lO"^® (9.378 x lO"®) ' 

C 2 - 3.18 

C 3 - 2.97 MPa-m^^^ (2.7 ksl-in^^^) > 
C, - 0.8 

Cj - 77 MPa-m^^^ (70 ksi-in^^^) > 


( 6 ) 
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dc/dM 


1/2 

When SI unite ere used, AK end K ere given In MPe-m end 

eii nuix 

le given In s/cyele. When U.S. Cuetoaery ere uccd, &K .. end K ere 

€Zt OUIX 

1/2 

given In kel->ln end dc/dM le given In In/ cycle. 

Figure 4 shows a comparison between experimental growth rates [l3] and 

growth rates calculated with equation (1) for the 2219-T851 aluminum alloy 

material. These constant- amplitude data were for R ratios between -1 and 

0.7, and maximum-stress-to-yleld-stress ratios of 0.15 to 0.77. Data generated 

when the net-section stress exceed the yield stress were not included. For 

the 2219-T851 aluminum alloy material, the yield stress (0 ) was 360 MPa 

ys 

(52 ksl) and the ultimate tdnsile strength (a^) was 455 MPa (66 ksl). The 

dashed lines show factor-of-2 bands about the perfect agreement line. Most 

of the data were well within these bands. 

Figure 5 shows a plot of AK against dc/dN for several R ratios 

for 2219-T851 aluminum alloy material to Illustrate the sigmoidal shape of 

equation (1). The experimental data were obtained from reference [l4]. Ibe 

curves were calculated from equation (1) . The R ■ -1 data were obtained from 

a small center-crack tension specimen (W « 76.2 mm) and the other data were 

obtained from small compact specimens (W > 50.8 mm). The crack-growth 

coefficients (Cj^, C 2 , aad C^) used to calculate the curves were 

Identical to thohe shown In equations (6). However, the coefficient 

1/2 

for the small compact specimens was 38.5 MPa-m and for the small center- 

1/2 

crack specimen was 55 MPa-m . The coefficients (C^) were calculated from 

1/2 

the Two-Parameter Fracture Criterion [l5] using Kj. ■ 550 MPa-m'*’' and m *= 1. 

These values of Kp and m were obtained from the final crack lengths and 

maximum stress levels used in the constant-amplitude tests from reference [13], 


10 



(••• appendix C). The calculated curvee and the experiaental data are In 
good agreeaent. 

APPLICATION OF THE CRACK-CLOSURE MODEL AND RATE EQUATION 

The analytical crack-cloeure model and cracfc-groirth program (FAST - 
Fatigue Crack Growth ^alysla of Structures) was applied to constant-amplitude 
and aircraft-spectrum loading on 2219-T851 aluminum alloy sheet material. 

Under constant-amplitude loading, the model was exercised under simulated 
plane-stress and plane-strain conditions. Some typical examples of crack- 
surface displacements at maximum load and of crack-surface contact stresses 
at minimum load are given. Calculated crack-opening stresses are shown as a 
function of crack length, stress ratio, stress level, and constraint factor (a). 
The calculated crack-growth lives under various constant-amplitude lo Lng are 
compared with experimental data. 

Because most crack-growth life is generated at small crack lengths and 
low stess levels (plastic-zone size is small compared to sheet thickness) , 
nearly plane-strain conditions (a * 2.3) were found to give a good correlation 
under constant-amplitude loading. These conditions were also used to predict 
crack-growth life under aircraft spectrum loading. The crack-opening stresses 
were calculated from the model as a function of crack length and load history, 
and the crack-growth rates were predicted from equation (1) . The predicted 
lives under various spectrum loadings are compared with experimental data In 
the following sections. 


Constant-Amplitude Loading 

Crack-surface displacements and contact r. tresses .- The normalized crack- 
surface displacements, EV/Sd, at maximum applied stress under constant-aqilltude 
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loading are 8ho«m as a function of coordinate location, x/d. In figure 6(a). 
The shaded region shows the plastic and residual deformations at maximum 
stress after simulated crack extension. The residual deformations along the 
crack surfaces are of nearly uniform length. The newly created crack-tip 
element comes Into contact after very slight unloading. The normalized contact 
stresses (o^/a^) at minimum applied stress ■ 0) are shown In figure 6(b). 

The elements are In contact all along the cracV surfaces, but only the elements 
near the crack tip carry significant loads. The elements near the crack tip 
yield in compression, whereas those near x * d are In tension. For x 
greater than d, the stresses (not shown) are In tension but rapidly drop 
to the applied stress, ^ Increases. 

Crack-opening stresses .- The crack-opening stresses are calculated from 
the contact stresses at minimum load as shown In appendix A. In figure 7, 
some typical crack-opening stresses are shown as functions of crack length. 
Results are shown for three stress levels used in the constant-amplitude 
tests (R » 0). At low stress levels the crack-opening stresses rapidly 
reach a "stabilized" level. However, at the highest stress level (276 MPa) 
the crack-opening stresses did not stabilize but did oscillate about a mean 
level for mos; f the simulated test. The sharp changes in at certain 

crack lengths were due to the element "lumping" procedure described In 
appendix A. Basically, this procedure lumps adjacent elements together to 
form a single element to keep the total number of elements to reasonable size 
and to save computer time. (The sharp changes in can be eliminated by 

using smaller elements.) These sharp changes In have less than 1 percent 

effect on the effective stress range. In the constant-ampxitude crack-growth 
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rate corrslfttioo (fig. 4), tli* •Cablllssd craek-<9«alaf sttMSM or tbo 
avcrago voluoo (otaoiin by tho d ath o d llaoo) voro oood. 

The varUtlon of atabUixod (or avorafo) cradt-opaotag otrasoM vlcb 
attaaa ratio* E* for • 0^/3 la ahowa la flgora 8. Tba curaas 

aho» bow aarlaa aa conditiooa cbsaga froa alaulatad plaaa atreaa 

(a ■ 1) to alaulatad plaoe atraln (a • 3). At aoy E ratio, tba S /S 

O SftZ 

values are lower for higher values of a. At E • 0* the crack-opening 

stress was about 0 47 for plane stress and was about 0.25 for 

plane strain. These results are In good agreeaent with flnlte-eleaent 

analyses under plane-stress £63 and plaitt-straln conditions (by the author). 

At negative E ratios, appears to be a linear function of the E 

ratio. At positive R ratios, rapidly approaches the dashed 

line. For E ratios greater than about 0.7, S/S equals S . /S , 

o sny wxti SAX 

and the values are also independent of the constraint factor. 

The influence of the ■axtaMa applied stress level on crack-opoilng 
stresses axe shown in figure 9. At the lower E ratios, the 
applied stress had a drsautic affect on crack-opening stresses, but at an 
R ratio of 0.7 the naxlaua applied stress had no influence. The Influence 
of the aaxl»Hi applied stress was aore pronounced under plane-stress con- 
ditions (dashed curves) than under plane-strain conditions (solid curves). 

Crack-growth calculations .- Table 1 shows the experiaental lives (N^) 
and the ratio of calculated to e]q>erlaental crack-growth lives (Np/N^) fron 
the initia? crack length c^ to the final crack length c^, for the constant- 
aaplltude tests on 2219-T8S1 aluninun alloy material. The closure model 
used a *■ 2.3. The ratio 8p/N^ ranged froa 0.6A to 2.09. The lowest and 
highest values of were obtained from the tests conducted at the highest 
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•trass The tsn vslus of vss 1.06 sad Chs stsadsrd dsvlstloa 

was 0.34 for sU tssts shovn in t^ls 1, sacspt ths low strsss leoel (55 Ws) 
test with R ■ -0.3. 

Ths two tssts iadlcatsd, by footaots (b) ia tsbls 1, wars not ussd in 
obtaining ths crack-growth constants (<!|^ sad C 2 )* Ths 55 MPa and R • -0.3 
crack-lsngth-agalnst-cyclss curvs did not agrss with two othsr tssts in 
rsfsrsQcs £l33 at ths saas strsss Isvsl and strsss ratio. Ths 276 MPa and 
nsgatlvs R ratio (-0.3) tsst was not Includsd bscauss, during ths siaulsLsd 
tsst» ths calculatsd valuss did not stabillss with Incrssaing crack 
Isngth. 

Spsctrua Loading 

Crack-surf acs dlsplaceaents and contact strsssss .- Ths noraallssd ccack- 
surfacs displacsasnts at asriaua applied stress, and ths noraalised contact 
stresses at atntnua applied stress under a typical spectrun loading, are 
shown as a function of coordinate location in figure 10. The dlsplacenent 
profile is quite siailar to that shown for constant-aaplitude loading (fig. 6(a)), 
but snail differences in the displacoKnt profile cause large differences in 
contact stresses, as shown in figure 10(b). The calculated crack-opening 
stresses reflect the influence of these irregular contact stresses on sub- 
sequeat crack growth. 

Crack-opening stresses .- The variation of crack-opening stress with crack 
length for a typical spectrum loading test is shown in figure 11. The half- 
length of the elox notch (c^) was 3.2 am. The specimen was cycled under 
constant-aaplitude loading (S • 69 MPa) at R ■ 0 until the crack grew 
to a crack half-length (c^) of 3.8 an. Next, a typical filter spectrum was 
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applied to the speciaen. The ■arlana stress ms shoot 183 H8s and the ifii-if 
stress was shoot -30 MPs. The psrtlcolsr spectros loads applied are gives Is 
reference [l3] for test M91. The calculated crack-opening stresses plotted 
In figure 11 rtuiw only a ssall fraction of the nmher of values conputed fron 
the nodel. The crack-opening stresses follcnr a very Irregular pattern while 
the cyclic loads are applied; even ao, tlwy tend to oscillate about a 
nean value. 

The use of an "equivalent" crack-opening stress concept would greatly 
reduce the cosputer tines required to cmplete a slmlated test. The use of 
an equivalent stress Is Justified, because the crack-opening stresses 
stabilise under constant-aaplltude loading at low to nedlun stress levels, 
and they tend to oscillate about a nean value under spectrum loading. The 
equation used to calculate an equivalent crack-opening stress, S^, was 



(7) 


where the suasMtlon was performed over the crack extension Increnent 

c. + 5 p to c. -I- 10 p . The maximum plastic-zone size, p , was 
i ^max i ^max '^max* 

calculated using the maximum stress in the spectrum. For extrenely high 
stress levels and low R ratios, where values do not stabilize, the 

simulated test specimen may fail before the equivalent crack-opening stress 
routine is activated. The dashed line in figure 11 shows the calculated 
equivalent crack-opening stress. The predicted crack-growth life using 
was 3.5 percent less than the predicted life using S^, but the computer time 
was only about one-half as large (2.6 to 5.6 minutes). 


15 



Cr«ck-growth predlctloo». - Th« crack-growth rat*, at each load cycle. 


was coaputed froe aquation (1), using the current values of S , S . , 

MMX in 

and S^. Equation (1) predicts retardation (or acceleration) if is 

larger (or sealler) than the crack-opening stress that would have been pro- 
duced under constant-aaplitude loading at and To deaonstrate 

how crack-growth rates were calculated under variable-aaplitude loading, an 
example is given. Figure 12 shows a typical variable-aaplitude load history. 
The growth rate was coaputed froa equation (1) using 


AK - AS /irc F 

eff eff. 

k 


( 8 ) 


where is the effective stress range on the kth cycle. The growth 

increment per cycle is 


Ac, 


(9) 


(Note that S and F are held constant during a small growth incr.£]Bent Ac*. 


See the section in appendix A on Crat.k Extension and Approximations.) On 


the first and second tensile load excursion, S , to S 

min, 

k 



AS 


eff 


k 



S 

o 


( 10 ) 


where k ■ 1 or 2, respectively. This equation was proposed by Fiber [2]. 
However, on the third tensile load excursion, greater than S^, 
therefore the effective stress range was assumed to be 
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( 11 ) 


• [«, 




(S . - S ) 

lOk o 




1_ 


^her^ k • 3 and is the power on the growth law. Thus, the growth 

iu;renent, AC 2 + ACj, is slightly larger than ACj^, if 

use of equation (11) was necessary because no crack-growth law, when expressed 
in t:e DS of a power function (C^ # 1), would sun to the correct growth increnent 
under variable-anq>litude loading. For Instance, If the load excursion S 
.0 extrenely small, then the sum of growth increments Ac 2 and 

^ 2 : Uould be equal to the growth increment If was less than 

S , then the growth increment Ac., should be equal to the growth increment 

o J 

Ac . liquation (11) accounts for these limiting behaviors. Equation (11) is 

apf lied only when is greater tlian acv. only when the current 

"k 

ma> linum applied stress is higher than the highest maximum stress occurring 
since a stress e* .urslon crossed On the fourth excursion, ^^eff 

again, computed from equation (10). Ihe effective stress range on the 5th 
and /t.i excursion were, again, coinputed from equation (11). But on the 6th, 
dth, ani 9th excusion. 


AS 


if 


-- S 


may 


min. 


( 12 ) 


where k =* 6, 8, or 9, respect ivoiv. Note that S v s . Equation (12) 
* max. maxe 

O D 

w. also proposed by *^,lber [i?]. 

The ASTW Ci4. 06.01 Round Robin involved five basic aircraf t-type load 
spectra :pp*ied to center-crack tonsiv>n specimens. Three of the spectra were 
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•ach applied at three different acale factors (aaae ahapa apactrun with dif- 
farant scaling of the stresses), and the other two spectra were each applied 
at two different scale factors. There were thirteen different spectra 
loading tests. *^0 test speciaens were precracked * 55 MPa, R >■ 0) fron 

an initial EDM notch (c^ •* 3.18 ■■) to the initial crack lengths given in 
table 2 before the spectrun loads were applied. See reference £l3] for more 
details. 

Figure 13 coapares predicted and experiaental crack-length-agalnst-cycles 
curves for a typical fighter spectrua. The speciaens were subjected to the 
saae spectrua, but with three dlfferen*' scale factors (0.2, 0.3, and 0.4). 

The predicted results (curves) are in good agreeaent with the experiaental 
data (syaibols). 

Figure 14 coapares experiaental and predicted lives for the thirteen 
spectrum loading tests (solid symbols) . The predicted results were obtained 
from the closure model with a > 2.3. The ratio of ranged from 

0.66 to 1.48. The mean value of Np/N^ was 0.98 and standard deviation was 
0.28. The results from the constant-anqplltude tests (open symbols) are also 
shown for comparison. 

Table 2 shows that the ratio lower for higher scale factor 

tests trlthln each spectrum. This behavior was attributed to the changing 
constraint at the crack tip. All predictions were made under nearly plane- 
strain conditions; however, the higher stresses are more likely to produce a 
plane-stress condition (plastic-zone size is large compared to thickness) , 
especially during a high overload, and, consequently, cause more retardation. 
These results suggest the need for a variable constaint factor which would be 
a function of the plastic-zone-to-thickness ratio. 
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Table 2 also presents the crack-growth predictions subwitted to the 
ASIM E24.06.01 Round Robin [l6]. These predictions were nsde using the 
analytical closure nodel but under nearly plane-stress conditions. In the 
original model, the bar elements were assumed to yield in cospression at 
instead of as is in the present nodel. This assumption strongly 

influenced the calculated crack-opening stresses. The crack-growth rate 
equation was also different from that used herein. The rate equation used in 
the original model is given in reference [9]. Crack-growth predictions 
under simulated plane-stress conditions gave higher crack-opening stresses 
than those under simulated plane-strain conditions (see figs. 8 and 9) and, 
consequently, retardation effects are much stronger. This may account for the 
large ratios obtained on test H88 and M89, where retardation effects 

dominated. 

Table 2 also presents crack-growth predictions using linear-cumulative 
damage (LCD), that is, no load- interaction effects. The predictions were 
obtained from Johnson [l7]. The predictions are quite good and indicate that 
retardation and acceleration effects nearly cancel each other for most of the 
spectra. Only test M88 and M89 had a sufficiently high overload such that 
retardation effects would dominate the life. Test M93 and M9A were dominated 
by crack-growth acceleration as is evident from the large ratios for 

linear-cumulative damage. 


CONCLUDING REMARKS 


An analytical fatigue crack-closure model was developed and used in a 
crack-growth analysis program (FAST) to predict crack growth under constant- 
amplitude and aircraft spectrum loading. The model was based on the Dugdale 


19 



model, but modified to leave pla8tically**de£ormed material in the vake of the 
ao ancing crack tip. The model was used to calculate the crack-opening 
stresses is a function of crack length and load history under simulated plane-* 
stress aul plane-strain conditions. 

A crack-growth rate equation was developed in terms of Elber*s effective 
stress-intensity factor range, threshold stress-intensity factor range, and 
fracture t'oughness, to give the "sigmoidal** shape commonly observed in fatigue 
crack-growth rate data plotted against stress-intensity factor range. The 
five crack-growth constants in this equation were determined from consta*«t- 
amplitude data on 2219-1851 aluminum alloy sheet material* The equation 
correlated the constant-amplitude data over a wide range of stress ratios 
and stress levels quite well. 

The analytical closure model, under nearly plane-stiain conditions, 
and the rate equation were used to predict crack growth under aircraft 
spectrum loading on the 2219-1851 aluminum alloy material* The model predicts 
the effects of load interaction, such as retardation and acceleration. The 
ratio of predicted-to-experimental crack-growth lives ranged from 0.66 to 1.48 
in thirteen spectrum load tests. However, many of the spectra were not 
suitable discriminators among load- interaction models because a linear- 
cumulative d*.mage model made equally good predictions. But these tests did 
verity whether the models could adequately balance retardation and acceleration 
effcjts. The only exceptions to this beliavior were tests M88 and M89 where 
retardation effects dominated, and tests M93 and M94 where acceleration effects 
dominated. The analytical crack-closure model and the proposed crack-growth 
law predicted crack growth in ail spectrum load tests quite well. 
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GOVERNING EQUATIONS FOR ANALYTICAL CRACR-CLOSDRE HODBL 

The following sections give the governing equations for the stress and 
deformation characteristics of the analytical crack-closure teodel. The 
sections include a discussion on crack-surface displacements » plastic-zone 
size, residual plastic deformations, contact stresses, crack-opening stresses, 
and crack extension. 

A schematic of the closure model at maximum and minimum applied stress 
is shown in figure 3. Because of synnetry, only one quarter of the plate was 
analyzed. A breakdown on the components of the model and the coordinate 
system used are shown in figure 15. The plate had a fictitious crack of half- 
length d and was subjected to a uniform stress S. The rigid-plastic bar 
element connected to point j was subjected to a compressive stress o^. This 
element is in contact when the length of the element is larger than the 

current crack-surface displacement . The stress Oj is applied to make 
Vj * ^ j ’ equations which govern the response of the complete system were 

obtained by requiring that compatibility be met between the elastic plate and 
all of the bar elements along the crack surface and plastic-zone boundary. 

The displacement at point i is 


S f(xp 


■ S °j 


(13) 


j-1 


where f(x^) and g(Xj^,x^) are influence functions given by 
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G(x^,Xj 


where 




- iiLz.nl 

E 






s«c 


ird 

W 


( 14 ) 


gCx^.Xj) » G(x^,Xj) >f G(-Xj,Xj) 


( 15 ) 




( 16 ) 



for 


k « 1 or 2, 


( 17 ) 
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•* Wj and b 2 ■ Xj + . Appendix B explains how these equations were 

developed. The bar element at point 1 (not shown In fig. 6) has a length 
(or residual deformation) of L^. The compateblllty equation (V^ ■ L^) 
is expressed as 
n 

g(x^,Xj) ■ S f(Xj^) - for 1 - 1 to n (18) 

j-1 

subject to various constraints. One type of constraint is caused by tensile or 
compressive yielding of the bar elements and the other is caused by element 
separation (V^ ^ L^) along the crack surface. The method used to solve 
equations (18) subject to the constraints is described in appendix D. 

Plastic-Zone Size and Approximations 
The plastic-zone size, p, for a crack in a finite-width specimen, was 
determined by requiring that the finiteness condition of Dugdale be satisfied. 
This condition states that the stress-intensity factor at the tip of the plastic 
zone is zero and is given by 

(K)g + (K)^ = 0 (19) 

^o 


where 


(K), 


max 


y'wd .ec ^ 


( 20 ) 
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alld 



( 21 ) 


Equations (20) and (21) are discussed In appendix B. Solving equation (19) 
for d, and noting that p ■ d - c, gives 



For a “ 1, equation (22) reduces to the expression derived by Rice [l8]. 

In the model, the plastic zone was arbitrarily divided Into ten graduated 
bar elements. The aspect ratios, 2w^/p, were 0.01, 0.01, 0.02, 0.04, 0.06, 

0.09, 0.12, 0.15, 0.2, and 0.3. The smallest elements were located near the 
crack tip (x * c) . Doubling the number of elements in the plastic zone had 
less than a 1 percent effect on claculated crack-opening stresses. At the 
maximum applied stress, the plastic-zone size was calculated from equation (22). 
The length, of the bar elements in the plastic zone were calculated 

from equation (13) as 

10 

‘•1 ■ 'l * ®max <“) 

i '1 
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where f(x^) end g(x^»x^) ere given bp equations (14) end (IS), reepectlvely. 
The bar elenents act aa rigid wedges. The plastic deforaatlon (L^) changes 
only when an eleoent yields in tension (o. > our ) or compression (o. < ) . 

) — 0 j “ o 

The constraint factor was assumed to apply only during tensile loading. 

Figure 3(a) shows a schematic of the model at the maximum applied stress. 

Here the crack surfaces are fully open and the residual deformations do not 
affect the crack-surface displacements. 

The division of the plastic zone into a number of finite elc<t 
would allow for the eventual use of a nonlinear stress-strain curve with 
kinematic hardening instead of the rigid-perfectly plastic assumptions used 
herein. This may give a more realistic crack-closure behavior. The models 
developed in references [?], [S], [lO], and [ll] are not well suited for 
incorporating a nonlinear stress-strain curve. 


Contact Stresses at Minimum Load 

When the analytical closure model is unloaded to the minimum load 
(fig. 3(b)), the bar elements in the plastic zone unloid until some of the 
elements near the crack tip yield in compression (a^ £ conq>resslve 

plastic zone, w, is shown as the shaded region at the crack tip. Depending 
upon the amount of closure and constraint factor, u varies from one-tenth 
to one-half of the plastic zone p. The original cyclic Dugdale model (no 
closure) [18] predicts that w equals 0.25 p at R = 0 and a ■ 1. Thus, 
crack closure greatly affects the size of the compressive plastic zone. The 
elements along the crack surfaces, which store the residual deformations, may 
come Into contact and carry compressive stresses. Some of these elements may 
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also yield in compression (a^ < -a ^) . The stresses in the plastic zone and 
the contact stresses along the crack surfaces are calculated from equations (18) 
with S For elements that yield in compression, is set equal to 

-0^ and the lengths of the elements are set equal to the final displacements 
at those points. That is, 


n 


j-1 


(24) 


For elements not in contact (L^ < V^), 0^ “ 0. 


Crack-Opening Stresses 

The applied stress level at which the crac) surfaces are fully open (no 
surface contact) denoted S^, was calculated from the contact stresses at 
^min* have no surface contact, the stress-intensity factor due to an 
applied stress increment (S^ - is set equal to the stress-intensity 

factor due to the contact stresses. Solving for gives 


min 


M-1 

s 

j.ii 


2a 


[sin Q, 


- sln'*^ Bj^] 


(25) 


where 
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for k ■ 1 or 2, 


( 26 ) 


and c is tbe current crack length nlnus Ac*. The incieaent Ac* is 
o 

the width of el 4ent n» and its sig^ficsnce is discussed in t^e next section. 
If • 0 for J ■ 11 to n " 1 at tbe slnisw applied stress, then the 
crack is already fully open, and cannot be detenined froai equation (25) . 
The stress 9^ at the crack tip changer fron conprcssion to tension when 
the applied stress level reaches S^. 

Crack Extension and Approxieations 

The analytical clot ire aodel provides no inforeation about the aaount 
ot craciL growth per cycle. Crack growth is slnulated by extending the crack an 
incremental value at the aonect of maximuD applied stress. The amount of 
crack extension, Ac*, was arbitrarily defined as 

Ac* * 0.05 p (27) 

max 

where is the plastic-zone size caused by the maximum applied stress 

occurring during the Ac* growth Increment. The nundier of load cycles, AN, 
required to grow the crack an increment Ac* was calculated from equation (1) 


and the cyclic road history. Typical values of Ac* ranged between 0.004 
to 0.1 on, depending upon the applied Cvress level and crack length. The 
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coefficient (O.OS) wee eclecced after e syateaetlc study wee aede on ttie effects 
of the else of Ac* on crack-K>penlng streseee and on conputer tine. Snaller 
values of Ac* Increase the conputer tine. The choice of the coefficient 
in this range did not adversely affect the crack opening stresses under 
constant-anplltude loading but sons sequence effects were lost under varlshle- 
aaplltude loading. 

The resulting conputer tines «rere considered reason^le. Under constant- 
anplitude loading, a siaulated fatigue-crack growth test, fron the initial 
crack length to failure, required 0.5 to 5 minutes on a CDC-6600 conputer. 

Under spectrun loading, a slnulated test to failure required fron 2 to 
15 otinuces. The conputer storage requirements for the analytical closure 
model and crack-growth program was low (less than 70K). The larger computer 
times occurred on tests chat were conducted at low stress levels and high R 
ratios. But these computer times could be drastically reduced If an '’equivalent" 
crack-opening stress were used In the latter stage of a simulated test. The 
use of an equivalent crack-opening stress is Justified because cr* ~-.k-opening 
stresses stabilize very quickly under low stress levels and high R ratios. 

The equivalent crack-opening stress is discussed in the section on 
"Application of the Crack-Closure Model and Rate Equation." 

The simulated crack extension Ac* creates a new bar element at the 
crack tip. The length of the new element is equal to the displacesmnt at the 
crack tip before the crack extends. 

The calculated crack-opening stresses and the boundary-correction factors 
were held constant while the crack was grown under cyclic loading (cycle-by- 
cycle) over the length Ac*. The number of load cycles, AN, required to grow 
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the crack an increoenc Ac* was calculated froa equation (1) and the cyclic 

load history. Hhen the aua of the crack growth Increnents (Ac) equalled or 

exceeded Ac*, the analytical closure nodel waa exercised. If AN reached 

300 cycles, the nodel was exercised trtiether or not Ac* was reached. This 

Units the nuhber of cycles that can be applied before the nodel is exercised. 

The increnent Ac* was set equal to sunnatlon of Ac's. Thus, Ac* was 

less than or equal to that conputed fron equation (27) and the nusber of 

cycles ranged fron X to 300. During the cyclic growth conputations, the 

cyclic stress history was nonitored to find Che lowest applied stress before 

(S , ) and after (S . ) the highest applied stress level (S ). The 

mliib nin^ nax^^ 

application of the analytical closure nodel consisted of: 

1. Applying o:lniaum stress crack length c. 

2. Applying maxi mum stress S at crack length c. 

nax^ 

3. Extending crack an increnent Ac*. 

4. Applying rainliaum stress at crack length c -t- Ac*. 

a 

5. Calculating new crack-opening stress fron equation (25). 

6. Continuing cyclic load history. 

7. Calculating new Ac* from equation (27). 

8. Repeating process when crack extension reaches new Ac* or AN 

reaches 300 cycles. 

Steps 1 through 4 require solving equations (18) for the c'ltact stresses 
and the stresses in the plastic zone. The overall procedure is repeated until 
the desired crack length is reached or the specimen fails. The failure 
criterion used to terminate the simulated test is discussed in appendix C. 
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To kftep the nueber o£ eleaente to a reasonable else (20 to 30), a 
"lunplng" procedure was used. The Iwi^lng procedure coablnea adjacent 
elcaenta (1 and 1 1) together to fora e single eleaent If 




c - a 


i-fl 


■f Ac* 


( 28 ) 


Thus, eleaenta near the crack tip are not as likely to be luaped together as 
those that are away iron the crack tip. In the luaplng procedure* the width 
of the luaped eleaent was the sua of the widths for the two adjacent eleaenta 
and the length was the weighted average of the two: 


L 


^1-H 

'‘ i ^ ’'i+1 


(29) 
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EQUATIONS FOR SIRESS-lBTENSm FACTORS AMD 
CRACK-SURFACE DISPLACEIOMTS 

The equations for stress-intensity factors and crack-surface displace- 
■ents for a crack in an infinite plate tiare obtained from the literature [19] . 
These equations are auaified herein for fini*‘e plates. Soae of the i^prozinate 
equations tor finite plates were verified with boundary-colloc.&tl'ti; 
analysri 


Stress-Intensity Factors for a Crack 
in an Infinite Plate 

The stress-intensity factors for the configurations shown in figure 2, 
with W equal to infinity, were obtained fron reference [19]. 

Reiaote uniform stress .- The stress- intensity factor for a crack in a: 
infinite plate is 


K„ - S /Hd 


(30) 


Partially-loaded crack .- The stress- Intensity factor for the con- 
figuration shown in figure 2(b) is 


K /fd 

00 71 


sin 


' (t) - (t)] 


(31) 
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Cr«ck-Surfac« Diaplaceatots 
for a Crack In an Infinite Plate 

The crack-surface dlsplaceeents for the configurations shown in figure 4, 
again, with W equal to infinity were obtained from Hestergaard stress 
functions given in reference [19]. The following sections give the dls- 
placeaent equations for renote unifora stress and the partially-loaded 
crack. 

Remote tinlfora stress .- Hie crack-opening dlsplaceaents for the con- 
figuration shown in figure 2(a) with U • » is 

V. - ^7^ (32) 

for |x] £ d, where n ** V for plane strain and n * 0 for plane stress. 

Partially-loaded crack .- The crack-opening dlsplaceaents for the con- 
figuration shown in figure 2(b) with W ■ is 

- V(x) + V(-x) (33) 


where 


V(x) 


2 (1 - n^) 

XE 



+ 




for |x| £ d. 


(34) 
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Flnlc»>Width CorMctlon* 

The equations given in the preceding sections for the stress-lntensity 
factors and crack-surface displaceesnts (eqs. (30) to (34)) were for a crack 
in an infinite plate. But these quantities are Influenced by the finite 
width of the plate. Because there are no closed fora solutions for finite- 
width specimens t approximate equations are developed herein. The stress- 
intensity factor for a crack in a finite-width plate is 

K - F (35) 

where F is the finite-width correction for the particular loading condition. 
Noting that the elastic crack-surface displacements in the region of the 
crack tip are directly related to the stress- intensity factor, it is proposed 
that the same correction factor be used for displacements: 

' ■ if) ' 


where V is the displacement for a crack in a finite-width plate subjected 
to the particular loading condition. Equation (36) gives very accurate 
crack-surface displacements in the region near the crack tip, as expected. 

Remote uniform stress .- The approximate boundary-correction factor for 
a central crack in a finite-width plate subjected co uniform stress [l9,2l] 
as shown in figure 2(a) is 
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The crack-surface displaceaents, given by equations (32) » (36) » and (37) 
were coopared with calculated dlsplacoKnts fron reference [20]* which used 
the boundary-collocation aethod, and the results were found to be within 
2 percent of each ocher for any value of x for 2d/W £ 0.7. The crack- 
surface displacements near the crack tip were very accurate, as expected. 

Partially-loaded crack .- The approximate boundary-correction factor for 
a central crack in a finite-width plate subjected to partial loading on the 
crack surface (fig. 2(b)) was obtained from the Infinite periodic array of 
cracks solution [l9] and was modified herein by replacing the term 
U/vd tan (ird/U) by sec (ird/W). It was shown in reference [2l] that the 
"secent" term gave more accurate results for stress intensity than the 
"tangent" term. The resulting boundary-correction factor was 



^sec 


F 


xd 

W 


(38) 
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where 



( 39 ) 


When « 0 and bj - d (unlformlyatrcased crack surface), equation (38) 
reduces to equation (37) , as expected. 
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FAILURE CRITERION 

The failure criterion used to tenninate the simulated fatigue-crack 

growth teats was the Two-Parameter Fracture Criterion (TPFC) [15]* The 

criterion involves two fracture parameters, Kp and m, that are determined 

from the fracture data. Because no fracture data were reported in 

reference [l3], the final crack lengths (recorded) and the maximum stress 

levels in the fatigue-crack growth tests (constant-asiplitude, variable- 

amplitude and spectrum loading) were used to estimate the fracture parameters, 

Kp and m. Figure 16 shows a plot of the elastic stress-intensity factors 

at failure, against the maxlmvim applied fatigue stress level. At a 

given stress level, the large amount of variation was attributed to the final 

crack lengths reported being less than the critical values and the maximum 

stress levels in the variable-amplitude and spectrum loading tests were 

larger than the actual failure stresses. The fracture parameter, m, was 

taken as unity because fracture data on 2219-T851 aluminum alloy in 

reference [l5] showed that m was unity for specimen thicknesses less than 

25 mm. The specimen thickness In reference [13] was 6.35 mm. The elastlc- 

1/2 

plastic fracture toughness, Kp, was chosen to be 550 MPa-m to fit the 
upper bound for most of the test data. The solid curve In figure 16 shows 
the calculated values at failure from the TPFC for the center-crack 

tension specimens (W * 152.4 mm) used In reference [13]. These results 
lndlc..te that specimens tested at low stress levels fall at lower stress- 
intensity factors than those tested at higher stress levels. 
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The maximum value of from the aolid curve wee about 77 MPa-m' 

Thia la the value of the coefficient Cj uaed In aquation (1). 


1/2 
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GAUSS-SEIDEL ITERATIVE METHOD WITH CONSTRAINTS 
The method used to solve the linear system of equations (eqs. (18)) 
is the Gauss-Seidel Iterative method [22] with constraints added. This 
nmthod was chosen because the nunber of equations to be solved were variable 
(10 to 40). The constraints are of two types. One type is caused by 
yielding (tensile or compressive) of the bar elements; the other is caused 
by element separation along the crack surface. The equations, which govern 
the deformations along the crack surface and in the plastic rone, are 


n 

2 Oj g(x^,Xj) - S f(x^) - (40) 

J-1 


where S Is the current applied stress, is the residual plastic deforma- 
tion for element 1, n Is the number of elements, and are the contact 

stresses along the crack surfaces and the stresses in the plastic zone. 
Equations (40) are solved with the following constraints: 


For elements in the plastic-zone region 


(Xj > c). 


if 0 , > oio set o. * cxa 

Jo Jo 


(41) 
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if Oj < -a^ set - -0^ (42) 

These constraints are for tensile aud compressive yielding, respectively. For 
elements along the crack surfaces (x^ < c), 


if 


> 0 et Oj ■ 0 


(43) 


if a. < -a set - -o (44) 

j o j o 

These constraints are for element separation and compressive yielding, 
respectively. 

First, rewriting equations (4) in the form 


n 


s 

1-1 




'i- 


and solving for 0^ gives 


(45) 
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n 

0^ * (S '• 

J-l 

for 1 ■ l,2,...,n. The matrix, is diagonally dominant, that is, 

each diagonal coefficient g^^ is larger, in absolute value, than the 
magnitudes cf other entries in row 1 and column 1. Next, make Initial 
guesses tor and Insert these into the right-hand side of equations (46) . 

If the initial guesses tor are zero, the resulting new approximations 

are 


(a,), 

(Oj)i 


ts 

[S f 2 “ 1>2 “ ® 21^^®22 


[S 


f - 
n 


(a.)i 


j-1 





(47) 


The additional subscript on denotes the iteration number. The stresses 

(o^)j^ are checkec against the constraint equations (41) through (44), as 
soon as they are computed, and are updated, if necessary. Note that the 
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Mvly-coaputed valuM of «re alway* used in the ri^t~hand eides as 
soon as they are obtalnud. hesubstitutlng into the right'-hand 

Bides of equations (46), subject to the constraint equations, glues for 

the next iteration. This process is repeated until the changes in axe 
less chan sons snail preset value (0.02 O^). The recurrence equation is 


(ai)i 


i-1 


Sfi- 1-1-2 

J-1 


n 


-s 

j-i+1 


^®j^I-l *1J 


/g 


11 


(48) 


where I Is Che current IteraCicn nuaber. Typically, Che number of elements 
(n) ranged from 10 to 35 and the iterations ranged from 2 to 13. 
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TABLE 1.- COMPARISON OF EXPERIMENTAL AND PREOIC1SBO CRACK 
GROWTH LIVES UNDER CONSTANT-AMPLITUDE LOADING 


max' 

MPa 


0.01 
O.Dl 
O.OI 
0.2 * 

0.3 


-0.3 

- 0.1 

0.01 

0.3 


4.89 

3.87 

3.91 

3.55 

3.87 

4.13 

3.85 

3.87 

4.00 

6.67 

4.57 

4.06 

3.94 

4.06 

3.94 

3.83 

3.94 

3.94 


24.6 

17.6 

18.7 
52.4 

6.9 

45.7 

43.0 

50.0 

42.0 

38.3 

48.0 
43.6 

12.1 

21.9 

21.4 

14.9 

12.9 

23.5 


Test 

«T 

Cycles 

Closure Model 
**P 

(o » 2.3) 
"t 

129,239 

1.42 

144,000 

1.7l(a,b) 

175,000 

1.49 

352,406 

0.98 

245,000 

0.84 

9,950 

1.04 

12,706 

0.96 

13,300 

0.98 

12,180 

1.03 

7,856 

0.93 

11,175 

1.19 

22,950 

0.80 

158,700 

0.87 

251 

2.09<‘» 

469 

1 

1.12 

846 

0.83 

1,693 

0.64 

14,870 

0.89 


This test is suspect based on other tests in reference 13. 
These tests were not used in obtaining crack-growth constants 












TABLE 2.- COMPARISON OF EXPERIMENTAL AND PREDICTED CRACK GROWTH 

LIVES UNDER SPECTRUM LOADING 







Test 

*^T 

Closure Model 

LCD^ 

SPECTRUM 

TYPE 

SCALE 

FACTOR 

SPECIMEN 



a « 2.3 
N« 

(a) 

N« 

[17 

Kp 


mm 

mm 

Cycles 

P 

P 








«T 

«T 


Air-to-Air 

Fighter 

0.2 

M81 

4.06 

13.0 

115,700 

1.20 

1.01 

. . 2 i 

0.3 

M82 

3.81 

35.4 

58,585 

0.71 

0.67 

o.'j: ! 

0.4 

M83 

3.81 

23.3 

18,612 

0.66 

0.64 

0. ri? 

Air-to-Ground 

Fighter 

0.2 

0.3 

M84 

MSS 

4.00 

3.66 

55.9 

44.1 

268,908 

95,642 

1.48 

0.96 

1.47 

0.89 

1.09 

0.81 

0.4 

M86 

3.87 

32.8 

36,397 

0.72 

0.65 

0.70 

Navigation 

0.3 

M88 

3.81 

45.8 

380,443 

1.47 

2.52 

0.73 

Fighter 

0.4 

M89 

3.81 

38.4 

164,738 

1.06 

1.47 

0.62 

Composite 

0.2 

M90 

3.87 

51.6 

218,151 

1.14 

1.17 

1.09 

Mission 

0.3 

M91 

3.81 

36.1 

65,627 

0.85 

0.78 

0.87 

Fighter 

0.4 

M92 

3.81 

29.5 

22,187 

0.74 

0.69 

0.88 

Composite 

Mission 

0.2 

M93 

6.35 

13.6 

1,354,024 

1.02 

0.76 

\ 

3.00 

0.3 

M94 

6.54 

9.5 

279,000 

0.76 

1.25 

1.81 

Transport 











\ 

Mean 

0.98 

1.12 

1.12 






Std Dev 

0.28 

0.50 

0.63 


(a ) 

Predictions submitted in ASTM E24.06.01 Round Robin. 

^^^Linear-cumulative damage (no retardation and no acceleration) predictions from 
reference 17. 








2219-T851 (Ref. 13) 
t = 6.35 nin 

-1 ^ R ^ 0.7 


0.15 S ^ 



FIG. 4 - Comparison of experimental and calculated crack-growth 
rates for 2219-T851 aluminum alloy. 





(b) Contact stresses at minimum load 


FIG. 6 - Crack-surface displacements and contact stresses under 
constant-amplitude loading (R = 0) . 








ed crack-opening stresses as a function of R ratio under 
-eunplitude loading for simulated plane-stress and plane- 














FIG. 13 - Comparison of experimental and predicted 
curves for spectrum loading. 
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FIG. 14 - Comparison of experimental (N>jt) and predicted (Np) cycles 
to failure for 2219-T851 aluminum alloy material under 
constant-amplitude and spectrum loading. 






